Proposal for a loophole-free Bell test using homodyne detection 
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We propose a feasible optical setup allowing for a loophole-free Bell test with efficient homodyne 
detection. A non-gaussian entangled state is generated from a two-mode squeezed vacuum by 
subtracting a single photon from each mode, using beamsplitters and standard low-efficiency single- 
photon detectors. A Bell violation exceeding 1% is achievable with 6 dB squeezed light and an 
homodyne efficiency around 95%. A detailed feasibility analysis, based upon the recent generation 
of single-mode non-gaussian states, confirms that this method opens a promising avenue towards a 
complete experimental Bell test. 
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In their seminal 1935 paper, Einstein, Podolsky, and 
Rosen (EPR) advocated that if "local realism" is taken 
for granted, then quantum theory is an incomplete de- 
scription of the physical world ill]. The EPR argument 
gained a renewed attention in 1964 when John Bell de- 
rived his famous inequalities, which must be satisfied 
within the framework of any local realistic theory 0. 
The violation of Bell inequalities (BI), predicted by quan- 
tum mechanics, has since then been observed in many 
experiments 0, 0, 0. ' disproving the concept of lo- 
cal realism. So far, however, all these tests suffered from 
"loopholes" allowing a local realistic explanation of the 
experimental observations by exploiting either the low 
detector efficiency or the time-like interval between 
the two detection events 0- In this Letter, we propose 
a realistic experimental scheme based on the conditional 
generation of non-gaussian entangled light states and bal- 
anced homodyning, which can circumvent these difficul- 
ties and allows for a "loophole-free" experimental Bell 
test. 

A test of Bell inequality violation typically involves two 
distant parties, Alice and Bob, who simultaneously carry 
out measurements on two systems prepared in an entan- 
gled quantum state. The measurement events (includ- 
ing the measurement choice) at Alice's and Bob's sites 
must be spacelike separated in order to avoid the locality 
loophole, that is, to rule out any possible communication 
between the two parties. Optical systems are therefore 
particularly suitable candidates for Bell tests, because 
entangled photon pairs can now be generated and dis- 
tributed over long distances 0,0- However, the currently 
available single-photon detectors suffer from too low effi- 
ciencies rj, opening the so-called detector-efficiency loop- 
hole, that is, the experimental data can be explained by 
local realistic theories wherein the detectors only click 
with probability rj. This loophole is present in all optical 
BI tests today. 

In contrast, very high detection efficiencies can be 
reached in optical systems using balanced homodyne de- 
tectors [llJ, |llj, |l^ , opening a very promising alternative 



to Bell tests based on single-photon detectors. Several 
theoretical works have shown that a violation of BI may 
indeed be observed with balanced homodyning provided 
that some very specific entangled light states can be pre- 
pared [ij, llJ, [l^ . Although a violation up to the max- 
imum theoretical limit can be achieved with homodyne 
detection 0, all the states required in unfor- 
tunately appear to be experimentally infeasible. On the 
other hand, the experimentally accessible p^lpRllsj l two- 
mode squeezed vacuum states are unsuitable for a Bell 
test because they are characterized by a positive-definite 
gaussian Wigner function, which provides an explicit hid- 
den variable model for homodyne measurements. So 
far, no feasible experimental scheme has been found that 
could be used to prepare non-gaussian states that exhibit 
a violation of BI with balanced homodyning. 

In this Letter, we show that a strikingly simple optical 
setup (see Fig. 1) can be used to conditionally gener- 
ate non-gaussian states that are suitable for this pur- 
pose. Our scheme requires a pulsed source of two-mode 
squeezed vacuum state, which can be expressed in the 
Fock basis as 



(1) 



where A = tanh(r) and r is the squeezing constant. The 
basic idea is to "degaussify^' this state by subtracting a 
photon from each mode [ijl, 12(1 l2lj . More precisely, we 
produce a non-gaussian entangled state with the use of 
two unbalanced beam splitters BS^ and BSb with inten- 
sity transmittance T, followed by two photon-counting 
detectors PD^ and PDs such as avalanche photodiodes 
(APD) . The successful state preparation is heralded by a 
click of both PD^ and PD^. We shall see that a detector 
efficiency as low as 77 = 10% suffices for the preparation of 
states exhibiting a violation of BI. Very recently, the gen- 
eration of pulsed single-mode non-gaussian states by pho- 
ton subtraction from squeezed states has been demon- 
strated experimentally by some of us ^ . We therefore 
can realistically envision the experimental realization of 
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FIG. 1: Proposed experimental setup. The source (controlled 
by Sophie) is based on a master laser beam, which is con- 
verted into second harmonic in a nonlinear crystal (SHG). 
After spectral filtering (F), the second harmonic beam pumps 
an optical parametric amplifier (OPA) which generates two- 
mode squeezed vacuum in modes A and B. Single photons are 
conditionally subtracted from modes A and B with the use of 
the beam splitters BSa and BS b and single- photon detectors 
PDa and PDs. Alice (Bob) measures a quadrature of mode 
A (B) using a balanced homodyne detector that consists of 
a balanced beam splitter BS3 (BS4) and a pair of highly- 
efficient photodiodes. The local oscillators LOa and LOb are 
extracted from the laser beam by means of two additional 
beam splitters BSi and BS2. 

a Bell test with non-gaussian two-mode states prepared 
along the same lines. 

The photon subtraction can easily be understood in 
the limit of high transmittance of BSa and BSb {T ^ 1), 
since then the most probable event leading to a click at 
each PD is when a single photon is reflected by each BS. 
In this limit, the conditionally generated state is very 
close to a pure state |2Qj 

00 
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(2) 

where aA,B are annihilation operators. This state is non- 
gaussian, and therefore suitable for a Bell test based on 
homodyning. 

In the experiment proposed here, Alice and Bob mea- 
sure the quadratures — cos9 + sinOp^ and = 
cos^a;^ -|- sin0p^, which have continuous spectra and 
satisfy [x-', p*^] = iSjk- In order to use the standard 
Bell-CHSH inequality |2J], we discretize these quadra- 
tures by postulating that the outcome is -1-1 when a; > 0, 
and —1 otherwise. Alice and Bob must then choose ran- 
domly and independently between two different measure- 
ments, corresponding to the choices of two angles 81,62 



and 4'ii4>2- Their measurements outcomes are thus de- 
scribed by four variables ai, 02, fei, 62, with values -1-1 or 

— 1. One defines the usual Bell parameter 

S ^ (aifei) + (0162) + (02^1) - (02^2), (3) 

where {ajbk) denotes the average over the subset of the 
experimental data when Alice measured aj and, simulta- 
neously, Bob measured bk ■ As we shall see, exploiting the 
fact that PDa and PDb can be viewed here as "event- 
ready" detectors Q , one can prove that all local-realistic 
models for Alice and Bob measurements must satisfy the 
Bell-CHSH inequality \S\ < 2 23]. 

In the formalism of "event-ready" detectors introduced 
by John Bell (2^, one should know, by some initiating 
event, when a measurable system has been produced. 
The main idea is to pre-select - rather than post-select 

- the relevant events. For that purpose, one considers 3 
partners, Alice and Bob who perform the measurements, 
and Sophie who controls the source, see Fig. 1. The entire 
data analysis must then be performed on a pulsed basis, 
with Sophie sending time-tagged light pulses (local oscil- 
lator and squeezed light) to Alice and Bob. In each exper- 
imental run, Sophie records whether her photodetectors 
clicked, while Alice and Bob carry out measurements of 
one of two randomly chosen quadratures. After register- 
ing a large number of events, the three partners discard 
all events not corresponding to an "event-ready" double- 
click registered by Sophie. The correlation coefficients 
(ajbk) are then evaluated from all remaining events, and 
plugged into the S parameter Q. 

In a local realistic approach, the light pulses in each 
time slot carry some random unknown parameters /i, 
which will ultimately determine the sign of x^ and x^ . 
The crucial feature is that the experiment should be 
designed in such a way that the measurements at Al- 
ice's and Bob's sites, as well as the detection of chcks 
from the conditioning detectors by Sophie, are all space- 
like separated [l^. Therefore, the probability distri- 
bution p(/i) must be independent of the measurement 
phases ^1,2 and 0i,2- The measured sign sa on Al- 
ice's side (resp. sb on Bob's side) only depends on fj. 
and 9 (resp. (f> on Bob's side). One can thus write 
(ajbk) ^ J dfi p{p) SAidj, fJ-) SB{4>k, IJ-), from which the 
derivation of the Bell-CHSH inequality is very standard 
[2^ . The only requirement needed in the above reason- 
ing is local realism, the absence of any other assumption 
being the condition for a "loophole- free" Bell test. 

Let us now evaluate the quantum mechanical value of 
the S parameter in the proposed experimental scheme. 
We only briefly outline the main steps of the calculation, 
the details will be presented elsewhere |2^. We model 
realistic photodetectors (which have a limited quantum 
efficiency rj < 1, and cannot discriminate between one 
and multiple photon detection events) as ideal detectors 
preceded with a "virtual" beam splitter of transmittance 
rj. Ideal detectors respond with two different outcomes, 
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either a no-click (projection onto the vacuum) or a chck 
(projection on all states with at least one photon). Simi- 
larly, a balanced homodyne detector with efficiency rysHD 
is modeled as a perfect homodyne detector preceded with 
a "virtual" beam splitter of transmittance 77bhd- At the 
output of the OPA, the modes A and B are prepared 
in a two-mode squeezed vacuum state, and the auxiliary 
modes C and D are in vacuum state. The Wigner func- 
tion of the state of modes ABCD is a gaussian function 
centered at the origin, 

Win,ABCD{r) = —^—^===exp [-r'^J^^r] . (4) 
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This state is fully characterized by the covariance 
matrix 7^^- — {riVj -f TjTi) — 2{ri){rj), where r = 
{x^ , . . . , ,p^) is a vector of quadrature compo- 
nents. The input covariance matrix is 7in = 7in,AB®/cr>j 
where / is the identity matrix and 7in,A_B denotes the co- 
variance matrix of the two- mode squeezed vacuum JQ. 
After combining modes A and C (B and D) on an un- 
balanced beam splitter BS^i (BSb) with transmittance 
T, the modes C and D are detected by the photon- 
counting detectors PD^i and PD^, while modes A and 
B are sent to the homodyne detectors. The covariance 
matrix 7out of the mixed gaussian state of modes ABCD 
after passing through the beam splitters BS^ and BSb 
and the four "virtual" beam splitters modeling the im- 
perfect detectors is related to 7in via a gaussian com- 
pletely positive map [13, 7out = Sr,SBsrinS'^g,S^ + G, 
where ^bs describes the mode coupling in BS^i and 
BSb, and the matrices Srj — ^/tjbwdIab © \/riIcD and 
G ~ (1 — ?7bhd)-^a_b® (l~fy)-^ci5 account for the imperfect 
detectors. 

The Wigner function of the state of modes A and B 
prepared by conditioning on observing clicks at PD^ and 
PDs can then be expressed as a linear combination of 
four gaussian functions, 
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(det7o,t)-i/2 ^-^^ q,[dei{T,T,.cD)] 



(5) 
-1/2 



where P 

stands for the probability of success, r = 
{x^,p"^,x^ ,p^), and we have defined qi — 1, 
<Z2 = 93 = ~2, (74 = 4. The various matrices appearing 
in Eq. jSJ are obtained from 7(^\, which can be split 
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into four smaller submatrices as 7„„4. — j 

i CD 

four submatrices are then used to define 
J - Tab - o-'^J^ho'^^^ where Ti^cD = Pcu, 
r2,CD = Pcu + /c ffi 0_D, r3,c_D = ^cD + Oc ® , and 

^4,CD = TcD + IcD- 

After discretization of the quadratures, the correlation 
coefficient E[9j, (j)k) = {a.jbk) can be expressed as 



sign(4 x^jPixt^ , xl )dxf^ dxl , (6) 
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FIG. 2: Violation of Bell-CHSH inequality with the 
conditionally-prepared non-gaussian state. (a) One- 
dimensional cut of the Wigner function of the conditionally 
generated non-gaussian two-mode state (with A — 0.5, T = 
0.95, and 77 = 30%) along the line = -x-*, p"^ ^ p'^ = 0. 
Notice the regions where W is negative, (b) Bell parameter 
S as a function of the squeezing A of the initial two-mode 
squeezed vacuum. The curves are plotted for perfect detec- 
tors (77 = ?7BHD = 100%) with T = 0.9 (solid line), T = 0.95 
(dashed line), and T = 0.99 (dot-dashed line). The open cir- 
cles mark the points where TA = 0.57. (c) Bell parameter S 
as a function of the efficiency r) of the single-photon detectors, 
for AT = 0.57, 77BHD = 100%, and the same three T's as in 
(b). (d) Bell parameter S as a function of the efficiency t^bhd 
of the balanced homodyne detectors, for AT = 0.57, 77 = 30%, 
and the same three T's as in (b). 



where P{xf.,x^^) is the joint probability distribution 
of the two commuting quadratures xf, and x^^ , which 
can be determined as a marginal distribution from the 
Wigner function given by Eq. jSJl- The correlation co- 
efficient E (and therefore S) can then be expressed in 
a closed form by analytically integrating the resulting 
gaussian functions. 

The results plotted in Fig. 2 were obtained for the op- 
timal choice of angles 61 = 0, O2 = 7r/2, (pi = —tt/4, 
(f)2 — tt/4. Figure 2(a) shows that the Wigner function 
W is negative in some regions of the phase space, which 
is a necessary condition for the observation of a violation 
of BI with homodyne detectors. Fig. 2(b) illustrates that 
the Bell inequality |5| < 2 can indeed be violated with 
the proposed set-up, and shows that there is an optimal 
squeezing Aopt which maximizes S. A simplified calcu- 
lation assuming perfect detectors with single-photon res- 
olution [2^ predicts that XoptT w 0.57. It follows from 
Fig. 2(c) that the Bell factor S depends only very weakly 
on the efficiency r] of the single-photon detectors, so the 
Bell inequality can be violated even ii ij < 10%. In con- 
trast, the efficiency of the homodyne detectors ryBHD must 
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be above ~ 90%, see Fig. 2(d). 

The maximum Bell factor achievable with our scheme 
is about <Smax ~ 2.046 which represents a BI violation of 
2.3%. This is a small, but statistically significant viola- 
tion, which should be possible to observe experimentally. 
The required degree of squeezing to get close to Smax cor- 
responds to A « 0.57, i.e., approximately 5.6 dB, a value 
that has already been achieved experimentally 0, . 
Another important parameter is the transmittance T of 
the beam splitters BS^ and BSb, which must be as high 
as possible for maximizing S. However, the probabil- 
ity of successful state preparation can be estimated as 
P K, 77^(1 — T)^, so that it quickly drops when T ap- 
proaches unity. Thus there is a tradeoff between S and 
P which needs to be optimized by taking into account 
the statistical uncertainties of the data. 

In order to be more specific, let us consider the single- 
mode version of our scheme, which has already been im- 
plemented experimentally |22| . It is based on a com- 
mercial cavity-dumped titanium-sapphire laser, deliver- 
ing nearly Fourier-limited pulses at 850 nm, with a dura- 
tion of 150 fs and a repetition rate of 790 kHz. Squeezed 
vacuum pulses, generated by parametric deamplification, 
are sent through a beam splitter, and the reflected beam 
is detected by a silicon APD. Conditional on observing 
a click, the transmitted pulse is prepared in a single- 
mode non-gaussian state, which is measured by homo- 
dyne detector with overall efficiency ?7bhd ~ 75%. This 
experiment gives us useful estimates for a possible Bell 
test. First, the delay between pulses (1.2 /is) allows am- 
ple time for individual pulse analysis. A fast random 
choice of the analyzed quadratures can be performed us- 
ing electro-optical modulators on the LO beams, trig- 
gered for instance by digitizing the shot-noise of locally 
generated auxiliary beams. Switching times around 100 
ns, associated with propagation distances of a few tens 
of meters, seem quite feasible. The APDs can be trig- 
gered only when a pulse is expected, reducing the effect 
of dark counts to a negligible value. The intrinsic APD 
efficiency is about 50%, but the filtering used to select a 
single mode currently reduces the overall rj to less than 
5%, which should be improved for accumulating enough 
statistics. 

We can thus define a set of realistic parameter values 
that should be reached in a loophole-free Bell test : with 
77 = 30%, T = 95%, and A = 0.6, BI are violated by 
about 1% if the homodyne efficiency ryBHD is larger than 
95% (see Fig. 2(d)). With a repetition rate of 1 MHz and 
P « 2.6 X lO^'*, the number of data samples would be sev- 
eral hundreds per second, so that the required statistics 
to see a violation in the percent range could be obtained 
in a reasonable time (less than one hour). In addition, 
the electronic noise of the homodyne detectors should 
be 15-20 dB below shot noise, attainable with low- noise 
charge amplifiers. All these numbers have already been 
reached separately in various experiments, but attaining 



them simultaneously certainly represents a serious chal- 
lenge. Nevertheless, taking into account many possible 
experimental improvements, the existence of an experi- 
mental window for a loophole-free test of Bell inequalities 
can be considered as highly plausible. As a conclusion, 
it appears that, with quantum continuous variables, a 
reasonable compromise can be found between the exper- 
imental constraints and the very stringent requirements 
imposed by a loophole- free test of Bell inequalities. 
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